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Black carbon measurements 

Continuous flow analysis was based on a steady sample flow and in-line detection of BC and other 
chemical substances as described in McConnell et al. (2007). In the cold room, previously cut one meter 
ice core sticks of 3x3cm, are melted continuously on a heated melter head specifically designed to 
eliminate contamination from the atmosphere or by the external parts of the ice. The melted ice from the 
most inner part of the ice stick is continuously pumped by a peristaltic pump and carried to a clean lab by 
Teflon lines. The recorded signal is continuous, integrating a sample volume of about 0.05 mL, for which 
the temporal resolution depends on the speed of melting, ice density and snow accumulation rate at the ice 



core drilling site. For annual accumulation derived from the WAIS and Law Dome ice cores, we assumed 
-3.1 cm water equivalent uncertainty in each year's accumulation from short scale spatial variability 
(glaciological noise) which was determined from several measurements of annual accumulation in 
multiple parallel ice cores notably from the WAIS Divide ice core site (Banta et al., 2008) and from South 
Pole site (McConnell et al., 1997; McConnell et al., 2000). 

Refractory black carbon (rBC) concentrations were determined using the same method as in (Bisiaux et 
al., 2011) and adapted to continuous flow measurements as described by (McConnell et al., 2007). The 
technique uses a single particle intracavity laser induced incandescence photometer (SP2, Droplet 
Measurement Technologies, Boulder, Colorado) coupled to an ultrasonic nebulizer/desolvation (CETAC 
UT5000) Flow Injection Analysis (FIA). All analyses, sample preparation etc, were performed in a class 
100 cleanroom using anti contamination “clean techniques”. The samples were not acidified. 

Calibration of the incandescence intensity was performed by introducing size-selected rBC nanoparticles 
(graphite particles of known density) into the SP2 using a differential mobility analyzer (DMA, 
performed by Droplet Measurement Technologies) (Slowik et al., 2007). Previous studies have found that 
the SP2 response is independent of rBC particle morphology (Cross et al., 2010), rBC coatings (Moteki 
and Kondo, 2007) and that the DMA calibration using graphitic rBC is applicable to rBC nanoparticles 
emitted by combustion (Schwarz et al., 2010). 

External calibration of the nebulizer/SP2 interface was performed using pure, commercial rBC 
dispersions (Tokai Carbon, Aquablack 162) with a similar surface chemistry to that of hydrophilic rBC 
aerosols found in the atmosphere (Moosmtiller et al., 2009). The instrumental response of these standards 
(0.03xl0 3 to 40xl0' 3 pg/L with 9 standards) was compared with a hydrophilic rBC dispersion produced 
by an inverted diffusion flame and ozonation (Kirchstetter and Novakov, 2007) and were found to give a 
comparable response. The commercial rBC was also used to prepare quality control standards, which 
were determined along with blanks every -2 hours of continuous melting. Blanks were found below the 
method detection limit (~ O.OlxlO" 3 pg/L, 3a). rBC concentrations were quantified by first calculating the 



individual rBC masses using the DMA incandescence calibration and then summing the rBC masses over 
three second intervals (3 sec). The FIA peaks in the resulting 3 sec mass time series were integrated and 
the peak areas was used to construct an external linear calibration, which was subsequently used to 
quantify the sample peak areas. Thus, interpolation below the detection limit is not performed as in 
Kaspari 2011. However, variability due to different calibration methods may not influence the relative 
changes observed in the records presented here and which are interpreted qualitatively. Indeed, ice core 
replicates analyzed at the beginning of each experimental day were found very reproducible with standard 
deviation lower than -10%. 

Spectral analysis 

Multitaper method (MTM), Figure 3 and SI-3: 

This method is independent of the spectral power, and even small amplitude oscillations may be 
considered as significant. The significance of each spectral peak is determined through an F-test (>0.90, 
>0.95, >0.99). In our cases, with monthly resolution over 150 years, we use a number of tapers of 3, and 
resolution = 2. For typical length instrumental climate records, this offers a good compromise between the 
required frequency resolution for resolving distinct climate signals (e.g., ENSO and decadal-scale 
variability) and the benefit of multiple spectral degrees of freedom (http://www.spectraworks.com/Help/). 
We used in particular the high-resolution multitaper spectrum which is the weighted sum of the K 
eigenspectra, with a red noise significance (robust noise), reshaped at 95% with linear fit. 
Blackman-Tukey method ( Figure SI-4 and SI- 5) 

This is the classical method for spectral analysis. In particular, we used it to determine the phase 
relationship between two data sets (between two sites or between rBC and Na), using two different filters 
to get the best spectral coherence for any given frequency. This coherence is a function of frequency with 
values between 0 and 1, and is a fraction of a common variance between two time series x and y through a 
linear relation. Coherence coefficients are given with 3 levels of confidence (high, medium and low). 
According to the software developer (Paillard, 1996), coherence is considered non-zero when coefficients 



reach values > 0.38. Phasing between the same periodicities was also determined through this method, 
and is given as a positive or negative value in radians as a function of delay between the two series. 


Variance estimation 

For variance calculation, we used this following equation: 


n 

where x is the average of the data set (in log, thus the geometric mean) and n is the size of the data set. To 
represent the period of low variance (Fig. SI-6) we choose to apply a 21 months running variance to the 
data set. 
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Figure SI-1: 1970-2000 monthly (dots) and annual smoothing (line) for Na at WAIS and Law Dome. 
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Figure SI-2 
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Figure SI-2: Distribution of rBC (left) and Na (right) data points at WAIS (top) and Law Dome (bottom) 
(black dots) and corresponding modeled log-normal fit, for 1850-2001 and 1970-2001. 


Figure SI-3 
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Figure SI-3: Spectrum obtained with Multi Taper Method for WAIS (a) and LD (b) monthly Na records for 1850-2001 period and confidence 

levels in red lines (90, 95, 99%). Numbers on graph indicate significant periodicities in years. 


Figure SI-4 
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Figure SI-4: Top. Average (thick) and max-min (thin) coherence coefficients for 1970-2001 period between rBC and Na records for WAIS(a) and 
LD(c). Non zero coherence is > 0.38 and shaded. Time resolution is 0.4yr (Blackman Tukey spectrum using Bartlett window, bandwidth 0.098). 
Bottom: Phasing between rBC and Na periodicities in radians for WAIS(b) and LD(d). Na leads for values >0. ENSO bands (~0.15-3 cycle/yr or ~3- 
7 years) and 1-yr bands are shown as dash lines. Black numbers indicate delay between rBC and Na in years at those corresponding 

periodicities. 


Figure SI-5 
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Figure SI-5: Average coherence (thick) and max/min (thin) between LD and WAIS periodicities for rBC (a) and Na (b) for the 1970-2001 period at 
a 0.4yr time resolution (Blackman Tukey spectrum using a Bartlett window, bandwidth 0.098, non-zero coherence >0.38). ENSO bands (~0.15-3 

cycle/yr or ~3-7 years) and 1-yr bands are shown as dash lines. 




Figure SI-6 
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Figure SI-6: moving variance (on log values) calculated for 21 months, from 1850 to 2000 at WAIS (a, c) and LD (b, d), for rBC (a, b) and Na (c, 
d). A period of low variance from 1955 to 1985 is common to the two sites, but only in rBC record. 


